Insulating behavior in metallic bilayer graphene: Interplay between density 

inhomogeneity and temperature 
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We investigate bilayer graphene transport in the presence of electron-hole puddles induced by long- 
range charged impurities in the environment. We explain the insulating behavior observed in the 
temperature dependent conductivity of low mobility bilayer graphene using an analytic statistical 
theory taking into account the non-mean-field nature of transport in the highly inhomogeneous 
density and potential landscape. We find that the puddles can induce, even far from the charge 
neutrality point, a coexisting metallic and insulating transport behavior due to the random local 
activation gap in the system. 



would typically manifest a transport activation gap sub- 
stantially smaller than the intrinsic spectral gap (i.e. the 
energy band gap) unless the band gap is much larger than 
the typical puddle-induced potential energy fluctuations. 

Our theory is based on a physically motivated idea: 
In the presence of large potential fluctuations V(r), the 
local Fermi level, /z(r) = Ep — V(r), would necessarily 
have large spatial fluctuations [particularly when Ep < s, 
where s = V rms is the root-mean-square fluctuations or 
the standard deviation in V(r)], leading to a complex 
temperature dependence of transport since both metallic 
and activated transport would be present due to ran- 
dom local gap. Below we carry out an analytical theory 
implementing this physical idea. We will see that this 
physical idea leads to the possible coexistence of metal- 
lic and activated transport, which explains the observed 
temperature dependence of BLG transport. 

We start by assuming that the disorder-induced po- 
tential energy fluctuations in the BLG is described by a 
distribution function P(V) which V — V(r) is the fluc- 
tuating potential energy at the point r = (x, y) in the 2D 
BLG plane. We approximate the probability P(V)dV 
of finding the local electronic potential energy within a 
range dV about V to be a Gaussian form, i.e., P(V) — 
, 1 | exp(— V 2 /2s 2 ), where s is the standard deviation 



Recent experiments [l|-[6( have revealed an intriguingly 
strong (and anomalous) "insulating" temperature depen- 
dence in the measured electrical conductivity of bilayer 
graphene (BLG) samples, not only at the charge neutral- 
ity point (CNP) where the electron-hole bands touch each 
other (with vanishing average carrier density), but also 
at carrier densities as high as 10 12 cm -2 or higher. ("In- 
sulating" temperature dependence of conductivity a(T) 
simply means an increasing a with increasing tempera- 
ture at a fixed gate voltage, which is, in general, consid- 
ered unusual in a nominally metallic system where the re- 
sistivity, not the conductivity, should increase with tem- 
perature.) Such an anomalous insulating temperature 
dependence of <j(T) is typically not observed in mono- 
layer graphene (MLG) away from the CNP although the 
gate voltage (or equivalently, the density) dependence of 
MLG and BLG conductivities are very similar with both 
manifesting linear-in-density conductivity away from the 
CNP and an approximatelya constant minimum conduc- 
tivity around the CNP 0,11] • 

In this Letter we theoretically establish that this 
anomalous insulating BLG u(T) behavior is likely to 
be caused by the much stronger BLG density inhomo- 
geneity [9j (compared with MLG) which gives rise to 
a qualitatively new type of temperature dependence in 
graphene transport, namely, the intriguing coexistence 
of both metallic and activated transport, hitherto not 
discussed in the literature. We therefore predict that the 
observed temperature dependence of BLG <r(T) arises 
from the same charged impurity induced puddles in the 
system which are responsible for the minimum conduc- 
tivity plateau at the CNP [l(|. We provide an analytic 
theory which appears to be in excellent qualitative agree- 
ment with the existing experimental results. One direct 
prediction of our theory, the suppression of the anoma- 
lous insulating temperature dependence in high mobility 
samples with lower disorder, seems to be consistent with 
experimental observations. As a direct corollary of our 
theory, we find, consistent with experimental observa- 
tion |lll4l3{ . that a gapped BLG (with the gap at the 
CNP induced, for example, by an external electric field) 



(or equivalently, the strength of the potential fluctua- 
tion). Then in the presence of electron- hole puddles the 
density of states (DOS) is reduced by the allowed electron 
region fraction and given by D(E) — J_ DoP(V)dV = 
Z?oerfc(— _E/v / 2s)/2, where erfc(a;) is the complementary 
error function and Do = g s g v m/ '(2ttH 2 ) is the DOS in a 
homogeneous system, where m is the band effective mass, 
g s = 2 and g v = 2 are the spin and valley degeneracies, 
respectively. We have Do — 2.8 x 10 10 cm~ 2 /meV with 
the effective mass m = 0.033m e (where m e is the bare 
electron mass). Note that the tail of the DOS is deter- 
mined by the potential fluctuation strength s. 

Since BLG is a gapless semiconductor the electron den- 
sity at finite temperature increases due to the direct ther- 
mal excitation from valence band to conduction band, 
and this thermal excitation is an important source of 
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temperature dependent transport. Thus, we first con- 
sider the temperature dependence of thermally excited 
electron density. The total electron density is given by 



D(E)- 



dE 



(1) 



where ft = \/ksT and Ep is the Fermi energy. When the 
Fermi energy is zero (or at CNP) all electrons are located 
in the band tail at T = and the electron density in the 
band tail is given by no = n e (Ep = 0) = Dos/ \/2tt. 
Note that the electron density in the band tail is lin- 
early proportional to the standard deviation s. At finite 
temperatures we find the asymptotic behavior of no{T). 
The low temperature {ksT/s -C 1) behavior of electron 
density at CNP becomes 



n e {T) = n 



7T 



(2) 



Thus, the electron density increases quadratically in low 
temperature limit. For homogeneous BLG with the con- 
stant DOS the electron density at finite temperatures 
is given by n e (T) = Do\n(2)ksT. The presence of the 
band tail suppresses the thermal excitation of electrons 
and gives rise to the quadratic behavior. However, at 
high temperature the density increases linearly with the 
same slope as in the homogeneous system, i.e., 



n(T) ~ Do 



I s 2 
\n{2)k B T+-- 



8 (k B Tf 



(3) 



In Fig. QJa) we show the temperature dependent electron 
density at CNP for different standard deviations. 

In the case of finite doping (or gate voltage), i.e., the 
Fermi level away from CNP, Ep ^ 0, the electron density 
of the homogeneous BLG for s = is given by 



noe(T) = DoE F 1+tln 1 



-l/i 



(4) 



where t = T/Tp and Tp = Ep/ks- At low tempera- 
tures (T <C Tp) the thermal excitation is exponentially 
suppressed due to the Fermi function, but at high tem- 
peratures (T> Tp) it increases linearly. In the presence 
of electron-hole puddles (s ^ 0) we have the electron den- 
sity at zero temperature for the inhomogeneous system: 



n e (0) = D E F 



-erfc | _ 
2 \V2S 



-1/2S 2 



2ir 



(5) 



where s = s/Ep. At low temperatures (T <C Tp) the 
asymptotic behavior of the electron density is given by 



n e (T) = rc e (0) + D E f 



,-l/2S 
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The leading order term is the same quadratic behavior as 
in undoped BLG (Ep = 0), but the coefficient is strongly 
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FIG. 1: (Color online) (a) The electron density at CNP as a 
function of temperature for different s. At T = the density 
is given by no = Dos/\^2tv. (b) The temperature dependent 
electron density at finite Ef for different s. For s/Ef 7^ the 
leading order behavior is quadratic while at a = the density 
is exponentially suppressed, (c) Total electron densities (solid 
lines) and hole densities (dashed lines) as a function of Ef for 
two different s — 30 meV and 70 meV. The linear line repre- 
sents the density difference n — n e — — DqEf, which lin- 
early depends on the Fermi energy. The densities at the band 
tails are given by n e {EF = 0) = n^(_E_F = 0) = Dqs/\/2tt. 



suppressed by fluctuation. In the case of s > Ef, the 
existence of electron-hole puddles gives rise to a notable 
quadratic behavior [see Fig.QJb)]. At high temperatures 
(T > Tp) we find 



n e (T) = n 0e (T) + 



D E f 



s 2 T F 



(1 



2 T 



(7) 



At CNP {Ep = 0) electrons and holes are equally 
occupied. As the Fermi energy increases, more elec- 
trons occupy increasingly larger proportion of space [see 
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Fig. 12c)]. For Ep ^> s nearly all space is allowed to the 
electrons, and the conductivity of the system approaches 
the characteristic of the homogeneous materials. In the 
presence of electron-hole puddles, there is a possible co- 
existence of metallic and thermally-activated transport. 
When electron puddles occupy more space than hole pud- 
dles, most electrons follow the continuous metallic paths 
extended throughout the system, but it is possible at fi- 
nite temperature that the thermally activated transport 
of electrons persists above the hole puddles. On the other 
hand, holes in hole puddles propagate freely, but when 
they meet electron puddles activated holes conduct over 
the electron puddles. Carrier transport in each puddle 
is characterized by propagation of weak scattering trans- 
port theory [9]. The activated carrier transport of pro- 
hibited regions, where the local potential energy is V 
less (greater) than Fermi energy for electrons (holes), is 
proportional to the Fermi factor. If a e and ah are the 
average conductivity of electron and hole puddles, re- 
spectively, then the activated conductivities are given by 
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(V) = <r e exp\p(E F -V)], (8a) 
(V) = a h exp[P(V - E F )}, (8b) 



where the density and temperature dependent average 
conductivities (a e and ah) are given within the Boltz- 
mann transport theory [9] by a e = n e e 2 (r) jra and 
&h = "-/ie 2 (r)/m, where n e and rih are average electron 
and hole densities, respectively, and (r) is the transport 
relaxation time which depends explicitly on the scatter- 
ing mechanism [9j. 

Now we denote the electron (hole) puddle as region 
'1' ('2'). In region 1 electrons are occupied more space 
than holes when Ep > 0. The fraction of the total area 
occupied by electrons with Fermi energy Ep is given by 
p = P(V)dV. Then the total conductivity of region 
1 can be calculated 

<ri = - [ EF (a e +a { h a) )P(V)dV, 

P J — OO 



a e + — e 
2p 



V2s V2 



^ • ( 9 ) 



At the same time the holes occupy the area with a frac- 
tion q = 1 — p and the total conductivity of region 2 
becomes 



a 2 



1 

3 Je, 



(ffW + a h )P{V)dV 
^^erfcf^ + ^']. ,10, 



The <7i and ai are distributed according to the binary 
distribution. The conductivity of binary system can be 
calculated by using the effective medium theory of con- 
ductance in mixtures The result for a 2D binary 
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FIG. 2: (Color online) cr t (T) at charge neutral point for dif- 
ferent s. Inset shows the thermally activated conductivity as 
a function temperature. 



mixture of components with conductivity a\ and a^ is 
given by [14 j 



*t = (p-2) 



(a 1 - <7 2 ) + J {a x - a 2 ) 2 



Aa\a 2 
(2p-l) 5 



(11) 

This result can be applied for all Fermi energy. For a 
large doping case, in which the hole puddles disappear, 
we have p = 1 and 02 = 0, then Eq. (fTTj) becomes a = 
<7i, i.e., the conductivity of electrons in the homogeneous 
system. 

We first consider the conductivity at CNP (Ep = 0). 
The conductivities in each region are given by 



a\ = a e 



&2 = cr h 



1 + JL e^ 2 /2 erfc(/3s/ V2) 
2p 

1 + — ^ V / 2 erfc(/3s/ yfi 
2qr) 



(12a) 
, (12b) 



where 77 = rih/n e is the ratio of the hole density to the 
electron density. Since the electrons and holes are equally 
populated we have p = q = 1/2 and a e = ah, then the to- 
tal conductivity becomes a t = \/a\ai — o\ ■ The asymp- 
totic behavior of the conductivity at low temperatures 
(ksT <C s) becomes 



a t(T) = a e 



2k B T 2 (k B Tf 



(13) 



The activated conductivity increases linearly with a slope 
^2/nkB/s as temperature increases. Because s is typ- 
ically smaller in higher mobility sample, the high mo- 
bility samples show stronger insulating behavior at low 
temperatures. The next order temperature correction to 
conductivity arises from the thermal excitation given in 
Eq. ^ which gives T 2 corrections. Thus in low temper- 
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FIG. 3: (Color online), (a) a t {T) for E F = 55 meV and for 
different s. (b) <Jt(T) for s = 50 meV and for several Ef = 18, 
36, 55, 78 meV, which correspond to the densities n = 0.5, 



ature limit the total conductivity at CNP is given by 



a t (T) = a(0) 



2 k B T 7T 2 (k B T 
s 6 



(14) 



At high temperatures (k B T 3> s) we have 



nk B T 2(k B T) 2 



(15) 



The total conductivity due to the activation behavior ap- 
proaches a limiting value and all temperature dependence 
comes from the thermal excitation through the change 
of carrier density given in Eq. ([3]). Thus at very high 
temperatures (T s/k B ) the BLG conductivity at the 
charge neutral point increases linearly with a universal 
slope ln(2) regardless of the sample quality. In Fig.[5]we 
show the calcuated temperature dependent conductivity 
at charge neutral point. 

At finite doping (Ef > 0) the temperature depen- 
dent conductivities are very complex because three en- 
ergies (Ep, s, and k B T) are competing. Especially when 



k B T <C s, regardless of Ep, we have the asymptotic be- 
havior of conductivities in region 1 and 2 from Eqs. (|9]) 
and (JTUJ), respectively, 



C2 = &h 



l + JLe-i/M 3 ./! I 

2p Mns/t-l/s 



1 + ^e 
2q V 



-l/2s" 



1 



7T s/t + I/S 



, (16a) 
(16b) 



where s = s/Ep and t = T/Tp. The leading order cor- 
rection is linear but the coefficient is exponentially sup- 
pressed by the term exp(— E F /2s 2 ). This fact indicates 
that in the high mobility sample with small s, the acti- 
vated conductivity is weakly temperature dependent ex- 
cept around CNP, i.e. Ep < s. Since the density increase 
by thermal excitation is also suppressed exponentially by 
the same factor [see Eq. ©] the dominant temperature 
dependent conductivity arises from the scattering time 
0. On the other hand, for a low mobility sample with a 
large s, the linear temperature dependence due to ther- 
mal activation can be observed even at high densities 
Ep > s. 

In Fig. [3] we show the total conductivities (a) for a fixed 
Ep and several s and (b) for a fixed s and several Ep. In 
total conductivity the activated insulating behavior com- 
petes with the metallic behavior due to the temperature 
dependent screening effect. When s is small the activated 
behavior is suppressed. As a result the total conductiv- 
ity manifests the metallic behavior However, for large 
s the activated temperature dependence behavior over- 
whelms the metallic temperature dependence, and the 
system shows insulating behavior. 

Finally, we discuss three important issues: (1) The 
same physics, of course, also applies to MLG graphene, 
but the quantitative effects of inhomogeneity (i.e. the 
puddles) are much weaker since simple estimates show 
that the dimensionless potential fluctuation strength s 
(= s/Ep) is much weaker in MLG than in BLG be- 
cause of the linear (MLG) versus constant (BLG) DOS 
in the two systems. In particular, s B lg/smlg ~ 32\m 
where n = n/10 10 , and therefore s_blg 3> smlg upto 
n = 10 13 cm~ 2 . Direct calculations [9( show that the 
self-consistent values of s tend to be much larger in BLG 
than in MLG for identical impurity disorder. In very low 
mobility MLG samples, where s is very large, the insulat- 
ing behavior of temperature dependent resistivity can be 
observed at high densities even in MLG sa mplesJUGil]. 
(2) We have neglected all quantum tunneling effects in 
our consideration because they are unimportant except 
at very low temperatures. In particular, Klein tunneling 
is strongly suppressed in strong disorder [ijj . (3) In the 
presence of a BLG gap (A s ), the situation becomes ex- 
tremely complicated since four distinct energy scales (s, 
Ep, k B T , A g ) compete, and any conceivable temperature 
dependence may arise depending on the relative values of 
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these four energy scales. It is, however, obvious that any 
experimental measurement of the activation gap (A ) in 
such an inhomogeneous situation will produce A a -C A g 
unless A g 3> s. The system is now dominated by a ran- 
dom local gap arising from the competition among s, A g , 
and Ep, and no simple activation picture would apply. 
This is precisely the experimental observations |llMl3j . 
Work supported by ONR-MURI and NRI-NSF-SWAN. 
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